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On the basis of the phase diagram of the uniform system we calculate the density profiles of a 
trapped polarized Fermi gas at zero temperature using the local density approximation. By varying 
the overall polarization and the interaction strength we analyze the appearance of a discontinuity 
in the profile, signalling a first order phase transition from a superfluid inner core to a normal outer 
shell. The local population imbalance between the two components and the size of the various 
regions of the cloud corresponding to different phases are also discussed. The calculated profiles are 
quantitatively compared with the ones recently measured by Shin et al., Phys. Rev. Lett. 101, 
070404 (2008). 

PACS numbers: 



I. INTRODUCTION 

The field of ultracold two-component Fermi gases with 
population imbalance is a very active area of research 
which in recent years has attracted a great deal of in- 
terest both experimentally and theoretically [l| . The ex- 
periments are carried out with clouds of atoms confined 
in harmonic traps and the analysis of the measured den- 
sity profiles is a key diagnostic tool to investigate issues 
such as shell structures and the phase transition from 
a superfluid to a normal state of the gas. Important 
achievements of this technique have been the observa- 
tion of phase separation between a superfluid core and 
a normal external shell driven by the degree of polar- 
ization 0, [1] , by the interaction strength tuned across a 
Feshbach resonance [1, Q and by the temperature of the 

gas iai. 

In the unitary limit, corresponding to a diverging scat- 
tering length between the spin-up and spin-down compo- 
nents of the fermionic mixture, the density profiles of the 
polarized gas have been investigated in details in Ref . @ . 
This theoretical study is based on an accurate determina- 
tion of the equation of state of the strongly-interacting 
normal g function of imbalance obtained using 

quantum Monte Carlo methods. It provides predictions 
for the shape of the profiles, for the density jump at the 
boundary of the superfluid core and for the critical polar- 
ization when the system turns fully normal, which are in 
excellent agreement with the experimental findings of the 
MIT group Q. In the Bose- Einstein condensate (BEC) 
limit of small and positive scattering lengths, the polar- 
ized gas is predicted to behave like a mixture of com- 
posite bosons (the bound dimer molecules) and fermions 
(the unpaired atoms) [H, [H, [l^. The density profiles 
in this regime have been investigated theoretically in 
Refs. [Tl, and experimentally in the recent study by 
the MIT group Q , where the Bose- Fermi mixture model 
is quantitatively tested at the level of mean-field theory, 
including also higher order corrections. At finite temper- 
ature, the density profiles of a trapped polarized Fermi 
gas have also been the object of theoretical investigations 
based on self-consistent approaches both at unitarity jlTj 



and in the BEC regime 

The phase diagram of the uniform gas at zero tempera- 
ture as a function of polarization and interaction strength 
has been calculated using quantum Monte Carlo (QMC) 
techniques in Ref. [l8|. This study provides a precise 
determination of the equation of state of four different 
phases competing for the ground state of the system: (a) 
the unpolarized and (b) polarized superfluid gas and (c) 
the fully and (d) partially polarized normal gas. The 
quantum phase transition from the normal to the super- 
fluid state is first order and is accompanied by a region 
of phase separation where the two phases coexist in equi- 
librium. Only in the deep BEC regime, where the Bose- 
Fermi mixture model applies, the transition from a fully 
polarized Fermi gas to a miscible mixture of few super- 
fluid bosons in a Fermi sea becomes second order. 

In this article we use the knowledge of the energy func- 
tionals of the uniform phases (a)-(d) as a function of po- 
larization and interaction strength and using the local 
density approximation we calculate the phase diagram 
of a trapped gas and the density profiles of both spin 
components. We determine the conditions for the ap- 
pearance of a superfluid unpolarized core and of a jump 
in the density profile signalling the occurrence of the first 
order quantum phase transition from superfluid to nor- 
mal. The calculated density profiles are compared with 
the experimental ones of Ref. JSI for different values of the 
interaction strength, from the unitary to the BEC limit, 
and for different values of the polarization. In Sec. |TT] 
we introduce the equations of state of the four uniform 
phases (a)-(d) and in Sec. Illll we discuss the equations of 
mechanical and chemical equilibrium when the harmonic 
external potential is present. In Sec. IIVI we present the 
results for the phase diagram, the density profiles and 
for the radii of the various shells present in the cloud. 
Conclusions are drawn in Sec. IVl 



II. UNIFORM PHASES 

In this Section we introduce the uniform phases consid- 
ered in Ref. [ll] and we discuss their energy densities as a 
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FIG. 1: (color online). Phase diagram for trapped atoms at 
zero temperature as a function of polarization and interaction 
strength. The regions correspond to different shell structures 
(see text). Inside the two (red) solid lines a jump in the 
density profile marks the superfiuid-normal first order phase 
transition. Above the (blue) dashed line the unpolarized SFo 
superfluid inner core is absent. On the right of the dotted line 
the polarized SFp superfluid phase is absent, while on the left 
the partially polarized Npp normal phase is absent. 



function of the concentration of spin-down particles (the 
minority species) and of the interaction strength. The 
natural unit of energy is provided by the Fermi energy of 
the spin-up majority component 
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where fc^t is the Fermi wave vector of the spin-up parti- 
cles fixed by their particle density n^. For trapped sys- 
tems, where the density varies with position, the above 
equation defines the local Fermi energy. 

(a) Unpolarized superfluid gas (SFq). In this phase the 
spin-up and spin-down particle densities are equal n| — 
ni = n/2. For positive values of the s-wave scattering 
length (a > 0) the energy density can be written in the 
form 



n /3 
fsFoW2) = -£6+ i-n^Epi 



2G(l/kF]a) 



(2) 



In the above equation Cb denotes the binding energy of the 
molecule (ef, = —Ti^ /ma? for a zero-range potential) and 
G(l/fc^-fa) is a dimensionless function of the interaction 
parameter rj ~ l/kp^a. A suitable parametrization of 
G{r]), which well reproduces the QMC results from the 
BEC to the unitary limit (see Ref. ^Ji]), is provided by 
the following formula 



G(ry) 



uq + OL\ arctan(a2?7) (r/ < 0.699) 



(3) 



where we defined the function 6(77) = 5(0.60/77)/(187r)[l-|- 
128(0.60/r;)3/2/(15V67r3/2)], while the constants ao-4 



have the following values: olq = 0.434, a\ = —0.292, 
a2 = 2.90, a-i = 0.0129 and 04 = -0.0100. The func- 
tion G as defined above is continuous at r/ = 0.699 with 
continuous first derivative. For kp^a <^ 1 it reduces to 
the function e(?7), which corresponds to the energy, in- 
cluding the Lee-Huang- Yang (LHY) beyond mean-field 
correction, of a gas of composite bosons interacting with 
the dimer-dimer scattering length add = 0.60a. 

(b) Polarized superfluid gas (SFp). This phase is char- 
acterized by a density np = of pairs and a density 
i^A ~ "-T^^-i of unpaired atoms, such that n — 2np-\-nA- 
The concentration of the minority atoms is denoted by 
y — n^ln^. For Xjkp^a > 0.5 an accurate parametriza- 
tion of the equation of state is provided by the following 
energy functional 
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The first term on the right hand side is the energy den- 
sity Q of an unpolarized superfiuid with density 2np. 
The other two terms correspond to the kinetic energy 
of a gas with density ua of unpaired fermions and to 
the interaction energy, treated at the level of mean field, 
between unpaired and paired atoms parametrized by 
the atom-dimer scattering length aad = 1.18a. With 
the above energy functional one recovers the Bose-Fermi 
mixture model in the deep BEC limit corresponding to 
l/ZcFTfl > 1- 

(c) Fully polarized normal gas (Npp). In this phase the 
density of spin-down particles vanishes (n^ = 0) and the 
energy density coincides with the one of an ideal Fermi 
gas 
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(d) Partially polarized normal gas (Npp). This phase is 
characterized by the concentration x = ni jn^ of the spin- 
down particles. For small concentrations (x <C 1) the 
dependence on x of the energy functional can be written 
in the form of the Landau-Pomeranchuk Hamiltonian of 
weakly interacting fermionic quasiparticles 0, [l^, [ll|, [T^ . 
[ll,[ll,[ll 
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where A is the binding energy of a single spin-down 
quasiparticle in the Fermi sea of spin-up particles, m* 
is its effective mass and F represents the interaction 
between quasiparticles. These quantities all depend on 
r\ = l/kp^a. We used the following parametrizations, 
which well reproduce the QMC results when 77 > 0: 
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where the first term on the right hand side is conveniently 
introduced in analogy with ^ and the constants P0-3 
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FIG. 2: (color online). Density (upper row) and column density (lower row) profiles at P = 0.8 for different values of 
the interaction strengtii. Column (a) corresponds to unitarity (l/fc^^a = 0), column (b) to l/fc^^a — 0.4, column (c) to 
l/k%^a = 0.8, column (d) to l/k%^a — 1.2 and column (e) to l/k%^a — 1.7. The solid (green) lines refer to the majority 
spin-up component, the dot-dashed (blue) lines to the minority spin-down component and the dashed (black) lines to the 
density difference — . The jump in the density n| of the minority component is shown with a vertical double arrow. The 
reference density no, column density rico and radius Ro all refer to a non-interacting Fermi gas with particles. 



have the following values: /3o = 0.986, /3i — 1.11, P2 — 
-2.23, /33 = 0.847; 



m*{r]) 



70 + liV + 72^7^ + 73?/^ + 74??^ , (8) 



where the constants 70-4 are given by 70 = 1.11, 71 = 
0.255, 72 = 0.337, 73 = 0.0579, 74 = -0.128; and finally 



(9) 



where the constants Sq_2 have the values: 6q — 0.279, 
5i = 0.438, S2 = 0.277. The equation of state (ED only 
holds for values of the interaction parameter rj < 0.73, 
since for larger values the normal gas becomes fully po- 
larized with X = (see Ref. [ill). We also notice that 
the values of the effective mass m* are smaller than the 
exact diagrammatic Monte Carlo results of Ref. [l^ (see 
comment in Ref. This inaccuracy is compensated 

in Eq. ^ by the coefficient F of the term, which fi- 
nally provides a very accurate fit to the QMC results at 
finite concentration. However, the reported value of F 
in Eq. ^ underestimates the interaction effects between 
quasiparticles. 



III. LOCAL DENSITY APPROXIMATION 

In the presence of the harmonic confinement 

Vhoir) = ^mu;lx^ + ^mujly^ + ^mujlz'^ , (10) 

the particle density depends on position and the concen- 
tration n^/rif varies locally. The most useful parameters 
for the phase diagram are expressed in this case in terms 



of the total number (-^x) of spin- up (spin-down) par- 
ticles. These are the overall polarization 



and the interaction parameter 

1 1 aho/a 
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(12) 



where Uho — \/TT'/rn{uJxUjyUJzy^'^ is the harmonic oscil- 
lator length and kp-, and 



1° = n^{0) are, respectively, 
the Fermi wave vector and the central trap density of a 
non-interacting Fermi gas with N-^ particles. 

In the local density approximation (LDA) the energy 
density of the inhomogeneous system is approximated 
at each spatial position by the one corresponding to a 
uniform gas at the local value of the density. This ap- 
proximation becomes exact in the limit of large particle 
numbers, when the density changes slowly on the typical 
length scale associated with the external potential. We 
construct the grand-canonical potential O at zero tem- 
perature as the integral over space of the ground-state 
energy density £, which depends on the local interaction 
parameter 77(r) = [n|(r)/n°]~^/'^/(/cp|a) and the local 
polarization P(r) — [n|(r) — n|(r)]/[n|(r) -I- n|(r)], that 



J rfr{£:[77(r),P(r)] + K(r)+nx(r)]T4„(r) 

-n|(r)^l -n|(r)^x} . (13) 

Here £ corresponds to the uniform phase present at posi- 
tion r, with the densities of the fermionic species replaced 



4 





1.0 


fit 


(a) 




0.8 






R 


0.6 








0.4 








0.2 


\ 






0.0 







l.Or 

0.8 



(b) 




0.4 0.6 
P 



0.0 
0.0 



0.4 0.6 
P 



1.0 

0.8 

R 0.6 

0.4 

0.2 

0.0 
0.0 



Rr 


(c) 











1.0 
0.8 
R 0.6 



0.2 



0.4 0.6 
P 



0.8 1.0 



0.0 
0.0 



Rr 


(d) 


^^ - — . — ___^SF=^i 













0.2 



0.4 0.6 
P 



0.8 1.0 



FIG. 3: (color online). Radii of the spin-up (R^) and spin- 
down ( ) components and of the unpolarized superfluid core 
Rc as well as of the superfluid shell Rsf as a function of 
polarization for different values of the interaction strength. 
Panel: (a) l/fc^^a = 0, (b) l/k°p^a = 0.5, (c) l/k%^a = 1, 
(d) l/fc^|a = 1.5. All the radii are divided by the radius Ro, 
which refers to a non-interacting Fermi gas with A''| particles. 



by the local densities, while /i| and are the chemical 
potentials of the spin-up and spin-down particles, intro- 
duced as Lagrange multipliers. When studying the su- 
perfluid phase we introduce the chemical potential fis 
of the pairs and the chemical potential fj,A of the excess 
spin-up fermions; chemical equilibrium at the interface 
between the superfluid region and the normal region is 
imposed by setting /15 = /i| -f /i| and /z^i = /i| . 

The density profiles are obtained by minimizing the 
grand-canonical potential with respect to the den- 
sities n|(r) and nx(r) of the two spin components. The 
condition of stationarity gives the LDA equations for the 
superfluid phase: 



^5 = dSsFp /dnp + 2Vho{r) , 
fj.A = dSsFp/duA + Vho{r) , 



(14) 



and for the normal phase: 



/i| = 9fNpp/9f^T + Vho{r) 
/ij, = dE^^^/dni + Vho{r) 



(15) 



Notice that mechanical equilibrium is enforced by the 
continuity requirement of the grand-canonical potential 
density along the cloud profile. For each point in the 
harmonic trap we solve the two set of coupled equa- 
tions, choosing the phase that minimizes the local grand- 
canonical potential and checking that the compressibility 
is positive. As usual, the integral of density over space 
allows us to fix the chemical potentials once we know the 
number of particles, or in our case, the total polarization 
and the interaction parameter. 



IV. RESULTS 

The shell structure of the atomic cloud for a given value 
of the polarization P and interaction strength 1/kp-^a is 
schematically described by the phase diagram shown in 
Fig. [TJ The various regions denote different sequences of 
phases starting from the center of the trap. The most ex- 
ternal shell is always composed by the majority spin-up 
atoms only (Npp phase). Outside the region delimited 
by the two (red) solid lines the sequence of phases cor- 
responds to a continuous density profile, while inside it 
corresponds to a density profile which exhibits a jump 
marking the first order phase transition from the super- 
fluid to the normal state. 

For large values of P and close to the regime of reso- 
nant interaction, the system is normal and displays an in- 
ternal shell of the Npp phase. The corresponding profile 
for both the spin-up and spin-down component is shown 
in column (a) of Fig. [21 For the largest values of the inter- 
action parameter corresponding to the leftmost 
region of Fig. [1] the central part of the cloud is either a 
polarized SFp superfluid (see column (e) in Fig. [2) or, for 
smaller polarizations, an unpolarized SFq superfluid sur- 
rounded by the SFp phase. In both cases the transition 
from the superfluid core to the Npp normal external shell 
is continuous and, within the picture of Bose-Fermi mix- 
tures, corresponds to a vanishing density of composite 
bosons. The superfluid-normal transition is instead first 
order in the case of the other three shell sequences shown 
in Fig. [TJ These are: the SFq ^Npp structure shown in 
column (b) of Fig. [1 the SFq-I-SFp ^Npp shown in col- 
umn (c) and the SFp — s-Npp shown in column (d). The 
jump occurring at the transition in the density profile of 
the minority spin component is shown in Fig. [2| with a 
vertical double arrow. 

The (blue) dashed line in Fig. [T] separates the region 
where an inner core of unpolarized SFq superfiuid is 
present (below the line), from the region where the polar- 
ized SFp phase survives at the center of the cloud. The 
dotted line, instead, separates the region where the first 
order phase transition is from an unpolarized superfluid 
to a partially polarized normal gas (right of the line), 
from the region where it is from a polarized superfluid to 
a fully polarized normal gas (left of the line). The values 
of P and 1/kp^a for which a first order phase transition 
directly between the SFp and the Npp phase is predicted, 
would correspond to a very tiny region around the dotted 
line in Fig. [TJ smaller than in the homogeneous case (see 
Fig. 4 of Ref. 18]) and involving such small densities of 
the minority component to be experimentally irrelevant. 

The structure of the cloud as a function of polarization 
and interaction strength can also be analyzed in terms of 
the radii of the various shells as shown in Fig. [3] The radii 
of the spin-up and spin-down components, respectively 
i?f and Ri , as well as the radii of the unpolarized super- 
fiuid core Rc and of the superfluid component Rsf are 
presented as a function of the polarization P for differ- 
ent values of the interaction strength. Panel (a) of Fig. [3| 
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FIG. 4: (color online). Comparison between calculated and experimental density (upper row) and column density (lower 
row) profiles for different values of polarization and interaction strength. Column (a) corresponds to P = 0.43 at unitarity 
il/k°p^a = 0), column (b) to P = 0.52 and l/k%^a = 0.35, column (c) to P = 0.57 and l/k%^a = 0.61, column (d) to P = 0.61 
and l/k'^ya — 1.2 and column (e) to P = 0.73 and l/k'^-^a = 1.6. The theoretical profiles are plotted as in Fig. [S] The 
corresponding experimental profiles are shown with thin lines. The reference density no, column density rico and radius Ro all 
refer to a non- interacting Fermi gas with N-f particles. The experimental profiles are from Refs. 



shows the behavior of the radii at unitarity (l/fc^^a = 0): 
the superfluid core exists here only if P < 0.77, while for 
larger values of P a partially polarized Npp gas is present 
for i? < At = 0.5 [panel (b)] the inner part 

of the cloud is always superfluid and a Npp normal shell 
still survives surrounding the SFq superfluid core. For 
larger values of [panels (c) and (d)] the super- 

fluid core is directly surrounded by the fully polarized 
Npp gas, but the radii Rsf and Rc are now different cor- 
responding to the structure of an unpolarized SFq inner 
core surrounded by a polarized SFp external shell. 

Finally in Fig.[3]we report the results of the compari- 
son between the calculated density profiles and the ones 
measured in the recent experiment of Refs. 0, 1^. We 
compare both density and column density profiles for dif- 
ferent values of polarization and interaction strength. It 
is worth stressing that there are no free parameters in 
this comparison, once the values of P and are 
extracted from the experimental data. The agreement 
is remarkable for all values of the interaction parameter 
from the unitary limit to the deep BEC regime, where it 
supports the picture of an interacting Bose-Fermi mix- 
ture 0. Small discrepancies, mainly visible as a broader 
density difference profile, are found close to the unitary 
limit [see panel (b) of Fig. 0] and might be attributed to 
temperature effects [20| . 



V. CONCLUSIONS 

We have investigated the density profiles of trapped 
polarized Fermi mixtures at zero temperature. The start- 



ing point is the reliable determination of the phase di- 
agram of uniform systems obtained in Ref. using 
QMC methods, whose implications for trapped systems 
are derived by applying the local density approximation. 
The comparison with the profiles measured in experi- 
ments shows a good agreement from the unitary limit 
to the BEC regime, providing new physical insights into 
the strongly interacting polarized normal and superfluid 
phases recently realized with ultracold atoms. 
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